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Abstract 

We describe a new type of three material microstructures which we call wheel assemblages, that 
correspond to extremal conductivity and extremal bulk modulus for a composite made of two materials 
and an ideal material. The exact lower bounds for effective conductivity and matching laminates was 
found in (Cherkaev, 2009) and for anisotropic composites, in (Cherkaev, Zhang, 2011). Here, we show 
different optimal structures that generalize the classical Hashin-Shtrikman coated spheres (circles) . They 
consist of circular inclusions which contain a solid central circle (hub) and radial spikes in a surrounding 
annulus, and (for larger volume fractions of the best material) an annulus filled with it. The same wheel 
assemblages are optimal for the pair of dual problems of minimal conductivity (resistivity) of a composite 
made from two materials and an ideal conductor (insulator), in the problem of maximal effective bulk 
modulus of elastic composites made from two linear elastic material and void, and the dual minimum 
problem. 

1 Introduction 

This paper introduces a new type of optimal isotropic structure for multimaterial composites: the wheels. 
They remain Hashin-Shtrikman coated spheres but consist of not only concentric annuli, but also an 
annulus with spikes: radial layers that alternate infinitely fast. We show optimality of wheel assemblages 
for the simplest structure of this sort: Three-phase two-dimensional composites of optimal conductivity, 
in which one of the phases is either an ideal conductor or a perfectly insulating material. 

Here, we briefly summarize the development of optimal structures. More detailed descriptions can 
be found in [3 1241 1§1 [TD] . In the pioneering work [TJ , Hashin and Shtrikman found the bounds and 
the matching structures for optimal isotropic two-component composites, and suggested bounds for 
multicomponent ones. Milton [23] showed that the Hashin-Shtikman bound is not exact everywhere 
(it tends to an incorrect limit when mi — >■ 0), but is exact when mi is large enough, mi > g-m- Nesi 
[22j suggested a new tighter bound for isotropic multicomponent structures, and Cherkaev [9] further 
improved it and found optimal laminates that realize the bounds. These latter two bounds coincide in 
the Ccise k-i — oo that is considered here. The method for obtaining bounds is based on the procedure 
suggested by Nesi [22] that combines the translation method [7] [24] and additional inequality constraints 
[3]. Cherkaev and Zhang [10] derive bounds for anisotropic composites made from two isotropic materials 
and void. 

The geometry of optimal isotropic two-material micro-geometries is intuitively clear: the stronger 
material "wraps" the weaker one so that the weak material forms a nucleus, and the strong one - an 
envelope. This principle is reflected in the first "coated spheres" geometry of optimal structures found 
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by Hashin and Shtrikman, and in following development: Lurie and Cherkaev |19l 121) suggested optimal 
high-rank laminates, Vigdergauz [55], Grabovsky and Kohn [T3] and recently Lui [TS] suggested special 
convex oval-shaped inclusions, Milton [21] introduced a general method of finding optimal shapes of 
inclusions, and recent paper by Benveniste and Milton [4] investigated "coated ellipsoids". For multima- 
terial case, the picture is much more diverse. Milton [23] introduced parallel coated spheres, Lurie and 
Cherkaev suggested multilayer coated circles [5D], Gibiansky and Sigmund [TT] suggested "bulk blocks", 
Albin and Cherkaev proved the optimality of three-material "haired spheres" [T]. All these structures 
admit separation of variables when effective properties are computed. However, study of these struc- 
tures has been overshadowed by the convenience of multiscale laminates and most recent results have 
been proven for laminates. Particularly, in [9] [10], it is shown that proper laminates are optimal for the 
considered here problem and its anisotropic generalization, respectively. 

The isotropic structures - wheels described here are direct generalizations of Hashin-Shtrikman coated 
spheres. They are intrinsically isotropic and simpler than isotropic laminates, and they nicely illustrate 
the qualitative properties of optimal structures. The topology of optimal structures depend on prescribed 
volume fractions of the materials being mixed. They undergo two phase transitions when the volume 
fraction of the best material decreases. The optimal microstructures are found by the combination of 
effective medium method and a procedure suggested in [2]. Similar "haired sphere" structures are studied 
in [T] and it is shown that they realize three material Hashin-Shtrikman bounds for a range of parameters. 

Acknowledgment The research is supported by the grant from DMS NSF. 



2 Equations and Notations 

Consider a composite: The materials ki occupy plane domains fi;, i — 1,2,3 C R2 whose union forms a 
unit cell fi. The areas rrii — of corresponding to prescribed volume fractions of each material 

are fixed: mi + 7712 + ma = 1, mi > 0, and there are no other constraints on fi^. The conductivity of the 
composite is described by a system of equations that relate fields e and currents j 

V X e = 0, or e = — Vu, V • j = in Q, j = he in Qi (1) 

where it £ Hi{Q,) is a scalar potential. The potential u and the normal current j ■ n are continuous at 
the boundaries dik between Q.i and Sl^, so the following conditions hold: 

T • (fii - efc) = 0, n • (ji - = But diu, (2) 

where r and n axe the tangent and normal to dik, respectively. 

In order to determine effective properties of an isotropic composite, we subject it to two homogeneous 
orthogonal external loadings eoa and eoi, of equal unit magnitude, 

lim ea{x) = eoa = i\] and lim e(,(a;) = eob = 

|x| — ^oc \ U y >oo 

and calculate the sum of energies. A pair of corresponding conductivity equations ([T]) (denoted with 
subindices a and b) differ only in boundary conditions. This pair can be conveniently viewed as a single 
problem for a vector potential U = (wa,Wb)"^ and 2x2 matrices E — {ea\eh) and J — {ja\jb) in a cell Q,: 

V • J = 0, E = -VU in n, J = kiE in lim E{x) = Eq ^ I. (3) 



E{x) is f2-periodic. 
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Energy In each material ki, the energy density Wi — Wia + Wn, of the pair is defined as a sum of two 
energies Wia and Wij,, caused by the loadings eoa and eob, respectively, 

= + Wib = ifcTr(£;£^) 

where Tr denotes the trace. 

Here, we assume that materials are ordered: 

< fci < A:2 < fcs = oo 

and fca is an ideal conductor. In fJa, field E is zero and the current is not defined. The energy of the 
ideal conductor is 

... r if £ = ... 

The energy Wo{Eo) of the whole periodicity cell has the form 

Wo{Eo) = inf ^ y / W,{E)dx, 

It defines the effective conductivity fc* through the formula 

Wo(E) = ^TT{E,,E'^) = k, (6) 

(because 'Tr{EQ Eq) = 2). Effective conductivity fc« depends on ki,mi, and on subdivision Qi. The 
geometrically independent lower bound of the efi'ective conductivity B is the solution of the problem 



I inf / W^{E)dx 



^ = ||o,,l^luo. ^1 -in^>./ W.iE)d.] (7) 
where mi m2 and ma are fixed fractions of materials, mi + m2 + ma — 1. 



3 Optimal bounds and optimal fields 
3.1 Exact Bounds and optimal fields 

Lower bound for effective conductivity The problem of exact bounds for an isotropic three 
material composite was solved in |22l[9] and the problem of bounds for an anisotropic composite in [lOj . 
It was found that effective conductivity fc* of a described two-dimensional isotropic composite is bounded 
from below by the bound B: 

fc» > B[m\,m2) for all microstructures (8) 

The bound has the form, see [9] 

(Bl if mil < mi < I, 
B2 if jTiii < mi < 77112, (9) 
B3 if < mi < mi2. 
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where 



B2 = k2 + 2 — (l^J^f. (11) 
mi 



and the threshold values mn and mi2 are 

mil = , (\/m2 — m2), mi2 = ^ (^/m2 — m2) . (13) 

fe2 + Kl /C2 

The bound is a continuously differentiable function of mi and m2. 

Optimal fields The bounds obtained correspond to the following fields in each phase [9] : The super- 
conducting phase naturally corresponds to zero field, e — in Qs, and the fields in fci and k2 are 

• If mi > mil (Hashin-Shtrikman bound), then 

TrE^^Hi, dct£;>Oinni, E^- — 3— Jfi in ^2, (14) 

ki ki + K2 



2ki ki + k2 



• If mi2 < mi < mil, (Intermediate bound), then 



TyE = 2- detE' = in fii, S = — ^/ in f72. (15) 

mi ym2 

• If mi < mi2 (Small mi bound), then 

TrE^^H2, detE = in Hi, Tr S = ^ 7^2 in ^2, (16) 

Kl K2 

fcl fc2 

These fields are derived from sufficient conditions for optimality, but they may or may not obey the 
differential equation and jump conditions. It is shown in [9] [10], however, that they are optimal: there 
exist multiscale laminates (see Figure [S] below) that has exactly the same effective properties as predicted 
by bounds and the fields that satisfy (|14|) . psp . H16p . In the present paper, we show different wheels 
structures that have the same effective conductivity, and the fields in them satisfy the same conditions. 

3.2 Comments on bounds derivation 

Here we briefly comment on the method used to derive bounds, and see [221 [9] [10] for derivations from 
the viewpoint of convexification. To derive the bounds, the problem ([7} is reformulated as a multiwell 
variational problem 

J ^ inf f F{E,j)dx, F(£,7)= min [W^,(£)+7,], W.{E) = ^\\Ef (17) 

E /q i — 1,2,3 Z 
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subject to 

(a) E = \7u, or \7 X E ^0, (b) I E dx ^ Eo (18) 



J > max min 

teR s(x):aQ = l, d{x):da=0 



it 

where 71 are Lagrange multipliers that account for constraints — mi. The nonconvex multiwell 

"energy" F of an optimal composite depends only on curlfree fields E, it is derived using an optimality 
condition which says that the larger conductivity k corresponds to smaller magnitude ||-E|| of E. 

Translation method [191 1211 [7J [21] replaces this variational problem with a minorant problem for 
E{x) £ ^2(^1) omitting the differential constraint (a) in (|18p but adding instead an integral constraint of 
quasiaffineness 

t / det (E) dx = t det / {E)dx, (19) 
Jn Jn 

see [251 I27| . that is satisfied for all gradient E = Vu but obviously not for all E. Here i is a real 
parameter (Lagrange multiplier). Energy W{E) and the determinant deti? are conveniently expressed 
through variables s and d 

s"" = {Ell + E22f + {E12 - E2if, d^ = {En - E22f + (£12 + £21)' 

as follows: 

WiE)^^{s' + d'), det{E) = ^{s'-d'). 
The minimization problem ([17p is estimated from below by the problem 

J iFtva.nsli^i^)'d{x),'y„t)dx -tdet{Eo) 

for any real t. We compute, fixing t 

^IransI = .^^i^g {^^''^'''^} (20) 

^transl ^ ^^^^^ + tdet{E) = i [{h + t)s^ + {h - t)d^] + 7i (21) 

where each well wf^^''^^^ is the translated energy of a material plus the cost 7;. The L2(n)-minimizers 5* 
and d are free from any differential constraints, and are subject to only integral constraints 

/ s{x)dx—l, / d{x)dx — 0. 
in Jn 

The convex envelope CF^j.^j^^i{s,d,t) with respect to s and d represents a lower bound of the composite 
energy for any value of t. The translation bound is obtained by maximization of CF^^j.g^j^gi{s, d,t) with 
respect of t. Notice that CF^.^^j^gi{s, d,t) — —00 if \t\ > ki because p^^^^^^ (|21[l is a quadratic saddle 
function. The constraint \t\ < ki must hold independently of the volume fraction mi. Because of it, 
the translation bound is nonexact for small values of nii and tends to a wrong asymptotic value when 
mi 0, see the discussion in [231 [9l [10] . 

The localized polyconvexity method used here for bounds refines the translation method by accounting 
for an additional pointwise inequality (Alessandrino, Nesi) [3] 

det(£) = - > 0, Vx £ n (22) 

that states that the Jacobian of two independent potentials does not change its sign in SI. To account 
for this constraint, we define the translated energy as +00 if the constraint is violated 

Fi{s,d,t) = <{ 2 j2 3 n ' * = 



F3{s,d,t) 



+00 if - < 

if s = d = 0, 
+00 otherwise 
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Notice, that the wells Fi are positive and grow quadratically with s, d for all real values of parameter t, 
but Fi and F2 become nonconvex functions of s, d when t > k\ and t > k2, respectively. 
The convex envelope CF of multiwell translated energy-cost function F 

F = min{Fi, F2, F3} (23) 

is nonnegative for all t. Therefore the constraint \t\ < ki is now lifted and the bound B is the maximum of 
F{t) with respect of t £ Ri. The bound is better than that obtained by the equivalent of the translation 
bound in this context. 

One can show that the bound is optimal when values of t belong to the interval [fci , fc2] . The translation 
bound corresponds to t = fci. When t > ki, the well Fi{s, d, t) becomes nonconvex. To address this case, 
we first find the convex envelope CFi{s,d,t) of well Fi 

( W/transl if ^2 _ ^2 > Q, < t < fci 
CFi{s,d,t)=l kis'^ + ji if - > 0, t > fci (24) 
[ +00 if - < 

The above bounds ([9]) and optimal fields are obtained by analyzing the finite-dimensional optimization 
problem of constructing CF. Three cases correspond to optimal values of t being topt = k\ (Bl), 
topt G (fci,fc2) (B2), topt = k2 (B3), respectively. 

Comment One can show that a point on convex envelope of piece-wise quadratic wells corresponds 
to the paraboloid with a coefficient equal to harmonic mean of the coefficients of translated energies, 
see bounds Bl and B3. The irrationality in the expression of B2 is due to optimization with respect to 
translation parameter t. The details of this calculation are shown in [9l [10] 



3.3 Effective medium theory 

We show that the bounds above are realizable by a special structures constructed by means of the effec- 
tive medium theory. Consider the Hashin-Shtrikman coated spheres scheme with anisotropic multilayer 
inclusions (wheels) of the following geometry. Assume that an infinite plane is filled with a material 
fc* and has a circular inclusion of unit radius. The central part of the inclusion of radius ro, hub, is 
filled by a material k2- It is surrounded by an annulus ro < r < 1 with an axisymmetric anisotropic 
material Kint{r) whose properties may vary with r but not with 9. Assume also that the eigendirections 
of Kintir) are r and 9. 

Because of axial symmetry of the inclusion, there is no need to consider two orthogonal external 
loadings: They correspond to the same but rotated solution. Assume also that the potential U , gradient 
field E and current J have the representation 

U ^u{r)cos{9), E = e{r)cos{9), J = j (r) cos{9) , 

that corresponds to a homogeneous applied field. 

In polar coordinates the conductivity equations have the following form: Ohm's law 

jr — Krrer Je = KegEg, (25) 

where the field e = Vu is represented as 

d 1 d 

Cr = —u eg^-—u (26) 

or r or 

and the equilibrium of currents (V ■ J = 0) 

'r^r^^^^^+^^9^^-' (2^) 
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The boundary value problem for u is 



Lu{u) = ( ^^^'■^ ~ ) U = r < Tincl 



du 
dr 



The energy of the assembly is 



w 



Kr 



at r 



du 
dr 



du 
dr 



u{0) = 



H —u 



r dr 



(28) 
(29) 

(30) 



Effective conductivity The effective conductivity of the assembly is computed using the Hashin- 
Shtrikman effective medium scheme, as follows: The homogeneous linear potential uqx = uorcosO is 
applied at infinity, 

lim u{u, 6) = uo r cos 9 (31) 



Outside the inclusion, where k = k,, potential u{r) satisfies Laplace equation and is equal to 

B 

U — uor ~\ 

r 



(32) 



where B is a perturbation caused by the inclusion. At the boundary of the inclusion, potential u and 
the normal current j„ = ^Tr'^ ^''^ continuous. 

If fc* is chosen so that i? = 0, the inclusion becomes "invisible" for an outside observer. In this case, 
fc« is called effective medium because the whole plane can be filled with identical inclusions at all scales, 
without changing the potential at the observer's position. 

The effective conductivity can be conveniently computed from the energy of the assembly and com- 
paring it with the energy of the "invisible inclusion" with undisturbed potential uq t. 



Wo 



1 . 



du 
dr 



2 2 
U 



r dr — ktUg 



Notice that fc« is independent of magnitude of u. Let us set uo = 1; then the energy of the inclusion is 
equal to fc» . 



4 Optimal Wheel Assemblages 

4.1 Wheel assemblages W{2, 13). Intermediate regime 

The Wheel W{2, 13) structure (nucleus from k2 then radial spikes from fci alternated with trapezoids 
from ks = oo). Figure [1] realizes the intermediate bound B2 ((Ojl when mi £ [mii,mi2]. The geometrical 
parameters of the structure are as follows: The central circle (r < ro = ^/rn2) is occupied by material fc2. 
The outer annulus y'm^ < r < 1 is filled with infinitesimal radial strips of constant width occupied by ki 
and trapezoids occupied by fcg. The total thickness of fci-strips is denoted bro, the total area 6ro(l — ro) 
of them is equal to the fixed fraction mi: 

mi = 6ro(l — ro) — by^m^il — ^mi) 

therefore 
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Figure 1: Cartoon of optimal Wheels W{2, 13). Intermediate case mi2 < nii < niu. White zones correspond 
to fci, grey - to ^2, and black to fcs = cxd 



b varies with volume fraction mi, 6 G [&ii, 612] where 

bll = &|mi=?7iii = -; ; — bl2 = 6|mi=mi2 = ~, — • 

Kl + fc2 K2 

The potential u{r) is continuous, we set it(0) — 0, u{l) = 1. The jump condition for the currents is 
not considered because Krr = 00 in the middle layer and the current is not defined. We compute 

(j^, r < ro ( k2l r <ro 

1 ro<r<l K{r) = i Kr = (x, Ke = ki-^ ro < r < 1 
r r > I y ktl r > 1 

where ro = ^/rn2. 
The energy is 

or, substituting the value of u(r), 

W ^ k2 + ki — ; — = A:2 + fci — = B2 

b mi 

This shows that W{2, 13) assemblage is optimal for the intermediate regime. 

Transitional points when mi = mn, the expressions for bound Bl and B2 in ^ coincide, and 
when mi = mi2, the expressions for B2 and B3 in ^ coincide. 

4.2 Wheel assemblage W{2, 13, 1): Large mi (Hashin-Shtrikman bounds) 

In the interval mi > mn, the optimal symmetric structure W{2, 13, 1) (Whole wheel. Figure [2]) can be 
constructed as follows: circles of the Wheel Structure W{2, 13) where (mi = mn) is enveloped by an 
annulus made from k\. When mi increases, the relative fractions in the central part stay constant, but 
the outer annulus becomes thicker. 
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Figure 2: Cartoon of optimal laminates [9] that realize the bounds for isotropic three-material composite 



To prove optimality of the stated structure, we notice (see [Zl[2|) that the effective conductivity fees 
of a coated circles structure satisfies the equation 

' ^ ' (33) 



+ ki 

where fci is an enveloping material, knuci is the core material, and c is the added fraction of fci. This 
implies that if a structure k^uci satisfies the Hashin-Shtrikman bound Bl, the envelope of this structure 
by fci also satisfies it by virtue of algebraic form of the bound and coated circles formulas, see [2] [7] 
One can check that the optimal W{2, 13) structure satisfies both bounds Bi and B2 if mi is equal to 
the critical value mi = mn. Therefore, W{2, 13, 1) structure (Figure [5} satisfies Bl bound for all larger 
values of mi. 

4.3 Wheel assembly W{2,123). Small mi 

In the interval mi < mi2, the optimal structures are wheel W{2, 123), see Figure[3l These structures are 



Figure 3: Cartoon of optimal wheel structures W(l, 123) that realize the bounds for isotropic three-material 
composite 

similar to the basic wheel W{2, 13) but the material fci in the spikes is replaced by a laminate of fci and 
fc2. At the point r — ro, the volume fractions are bi and 62, respectively. Materials fci and fc2 occupy the 
areas 

mi = 26iro(l - ro), m2 ^ rg + 2b2ro{l - ro) 
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and we find 

2ro(l-ro) 2ro(l - ro) 

The continuity of potential at r = ro and the sufficient optimahty conditions require the representation 
of the gradients at r = ro 

'^^^^(1) ^^^^i(o) ^- = ;^|(o) (35) 

Notice, that the continuity of current does not produce any constraints on fields because the current in 
superconductor phase is not defined. 

The optimal value ro of the radius of central patch filled with k2 is computed as 

/mi m2 
ro = fc2 ( — + — 

Ki fc2 



and corresponds to the energy equal to 



H/ = fe+(^ + Ji) (36) 



that agrees with the bound _B3. Condition 62 > is equivalent to condition mi < mi2 



4.4 Asymptotics 

li is interesting to find out how the three-material optimal structures degenerate into two-material ones. 
The known optimal two-material structures are Coated Circles VK(2, 1) when (ms — >■ 0, M^(3,2) when 
mi — >■ 0, and 14^(3, 1) when (m2 — >■ 0. Asymptotics of three material wheel when m2 ^ 1 and ma <^ 1 
are shown in Figure 2] upper line. The structures degenerate into coated spheres. In the first case, the 
radius of the central sphere (equal to y/m^) goes to zero, and in the second the thickness of the annulus 
goes to zero. If mi = 0, the three-material wheel structure degenerates into W^(2,23). The radius ro of 
the inner circle filled with ^2 equals m2, and the area fraction is m2. The remaining fraction m2(l — m2) 
of k2 is located in spikes, and fca in the trapezoidal areas between spikes. This structure is equivalent 
(it stored the same energy and has the same effective modulus fc, ) to coated circles M^(3, 2) where fca is 
placed in the center circle and k2 in the surrounding annulus, see Figure |4l lower line. The fields in these 
two structures are different but both satisfy sufficient conditions and the energy is the same. 

5 Dual problem. Two-material-and-void composites of min- 
imal resistivity 

Structures similar to those described above also deliver minimal resistance of a composite made from 
two conductors and an insulator. Using the representation for any divergence-free vector j 

j^RVv, J?=(_°^ (37) 

and accounting the identity R^R = I, we present the energy of a material in the form 

identical to 1301 Here, p = ^ is the resistance (a reciprocal to conductivity) of a material. Because 
energy form is identical, the bounds (|9]) are valid where ki must be replaced b pi, and it is assumed that 
Pi < P2 < p3 = 00. 
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Figure 4: Upper Left: Optimal wheels, m2 <^ 1, Upper Right: Optimal wheels, <C 1, Lower line: 
Limiting optimal wheel structure and equivalent coating circle structure: mi — > 0. 



Optimal Wheels structures of minimal resistivity The intermediate regime (bound B2, (|9j) 
is realized by the Wheel W{2, 13) described above. This time, the trapezoidal segments filled with pa are 
insulators, conductivity in the circumferential direction is zero and all current goes sequentially through 
the spikes and central circle. For small mi, bound B3 is realized by the Wheels structure W{2, 123). 

Compare the currents in wheel W{2, 123) in the primary problem with an ideal conductor in trape- 
zoidal domains and the same wheel in the dual problem with an ideal insulator in these domains. In the 
primary problem, the potential within each trapezoid is constant and the current in the spikes flows in 
circumferential directions across the layers, from one trapezoid to the next one. The current density is 
constant because layers have constant thickness. Another portion of currents flows through the central 
circle. The whole assembly is equivalent to a circuit: The layers in spikes are equivalent to sequentially 
jointed resistors, and the central circle is equivalent to a parallel resistor. In the dual problem, the 
trapezoids are filled with an ideal insulator and the circumferential current is zero. The current in spikes 
fiows radially along the layers and then fiows thought the central circle. The equivalent circuit consists of 
two parallel resistors that represent layers in the spikes, and a sequentially joined resistor that represents 
the central circle. These circuits are dual, because the field in parallel resistors is constant as the current 
is in the sequential conductors, and vise versa. 

Large values of mi correspond to the Hashin-Shtrikman bound and Wheels W{2, 13, 2). Notice that 
in each point of the same Wheel, the fields in the primary and dual problems are mutually orthogonal. 
In primary problem, the current in the spikes zone flows between the domain of superconductor across 
the layers, while in the dual one it fiows along the layers parallel to domains of insulator. In the central 
circle, their currents are also orthogonal. 

5.1 Laminates and Wheels 

Let us compare the optimal laminate structure found in [9] and |10) . shown in Figure [5] with the wheel 
structures found here, shown in Figure [6] The laminate structures [10] are anisotropic, and they depend 
on the degree of anisotropy and on volume fractions. Isotropic laminates are shown in the upper line. 
They also contain disconnected domains of k2 where the field E is proportional to a unit matrix, and 
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(fcijfca) or {ki,k2,kz) laminates that join these domains and direct currents between or around them. 
Larger amount of fci transforms the structure by placing the inclusions into a matrix from fci; these 
structures realize the Hashin-Shtrikman bound. The wheel assemblages have similar features, but are 
simpler and easier to calculate. 

Comparing the wheel assemblage with anisotropic laminates for bounds B2 and B3, having small 
volume fractions of mi (the lower left and central fields in Figure we notice that the wheel structures 
correspond to these laminates if they are curled so that the lower layer becomes a circle and the upper 
layers become spikes and trapezoidal domains. The curled structure gains the axial symmetry and keeps 
the property of {k\,k^) or (k\,k2,k->,) laminates that direct currents between or around domains of k2. 




L(13,2,13,2) Region Al 





L(13,2,I3) Regions 



L(13,2,13,l,l) Region Dl 





L(123,2) Region A2 



L(13,2) Region C 



L(13,2,l) Regions D2, E 



Figure 5: Cartoon of optimal laminates [5] that realize the bounds for isotropic three-material composite. 
They depend on the degree of anisotropy (vertical axis, top line corresponds to structures close to isotropy, 
bottom lime, to strongly anisotropic structures) and on the volume fraction mi of the material ki (the 
volume fraction increases from left to right) 



B3. 




B2. 




Bl. 




Figure 6: Cartoon of optimal Wheels structures that realize bound Bl, B2, B3. 



A natural generalization of a wheel for anisotropic optimal composites would be a structure assembled 
from confocal ellipses instead of concentric circles and orthogonal hyperbolic layers instead of radial layers. 
Similar structures (excluding those with hyperbolic layers) are described in [4] . 
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6 Elastic 2D structures from two materials and void with 
maximal bulk modulus 

The technique described above can be used to find tlie bound for the effective buik modulus of a composite 
made from two linear elastic materials and void (plane problem). The elastic energy density Wei is 
expressed through the stress tensor a as 

W=^[^ia^f+r^{a2f] (38) 

where ai is the normalized trace and (72 is a magnitude of normalized Deviator (trace-free part) of a 
2x2 symmetric tensor a: 



1 rr tJii + a-22 / (o-ii - 0-22)^ + 2o-^2 

"^ = 71 ^^^71— ' "^ = V 2 

Notice that Trtr^ = 2 (tJi + 0-2). Here it is convenient to use elastic moduli k and ri, the reciprocals of 
the plane bulk and shear moduli, respectively, and they are expressed through the Young's E modulus 
and Poisson coefficient v as 

1-1/ 1 + P 

''=^- 

The moduli are assumed to be ordered as follows 

Hi < K2 < H3 = 00, rji < ri2 < rjs = 00, (41) 
The Hooke's law corresponds to the Euler equation for energy (|38p minimization and has the form 

£11 = (k + ?7)f7ii + (k — 77)0-22 

£22 = (k + 77)0-22 + (k — 77)011 
£12 = 277O12 



The equilibrium equations 
allow for an Airy representation 



V • o = 0, = 0^ (42) 

a^R^{Vy)R, y = Vx. (43) 

Because the rotated tensor o is a gradient of a vector y, we may reduce the problem to the previous 
case and establish the bounds for k, in the same way it is done for two-material composites in [?]■ The 
lower bound of effective k modulus (that corresponds to the upper bound for effective bulk modulus) is 
obtained similarly to the conductivity case using relation (|19|l and (|22|l . 
Translated energies are 

^transl. , ^ _ 1 / (K + t)(oi)2 + (77- t) (02)^+7, if det(o) > . . 

" 2 \ +00 if det(o) < ^ ' 

When t £ [771,772], the energy Fi{s,d,t) becomes nonconvex, its convex envelope CFi{s,d,t) is 

r VF^'"'^^' if det(o) > 0, < t < 771 

CFr{s,d,t)=l m+nxi^aif +-yi if det(o) > 0, t > 771 (45) 
[ +00 if det(o) < 

In this case, the minimal energy corresponds to condition deto = in fii, corresponding to uniaxial 
stress in the first material. 
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To derive the bound, set the average stress field to be proportional to a unit tensor, ((to)j = 1, 
(0-0)2 ~ 0' ^^'-^ estimate the energy of a composite from below, obtaining bounds similar to [221 



K* > -B, B = max 

We compute: 

• If mi > mil, then topt ~ rji and B — Bl, 



Bl = -r/i + 



«:i + ??i H2 +t 



nil m2 



(46) 



«:i + m K.2 + rji 



(47) 



(Hashin-Shtrikman bound) 



• If mi2 < mi < mil, then topt = ^ 1 + ;i) ^ ^ _ 



(Intermediate bound) 



• If mi < mi2, then t, 



opt 



B2 = K2 4 

: ri2 and B = B3 

B3 = -r/2 + 



mi 

^ (1 — ym2)(«:i + rii) 
mi 



mi 



+ 



m2 



^1 +rii 1^2+ m 



(small mi bound) 
where the threshold values mn and mi2 of mi are: 



/m2(l — -v/m2)- 



/m2(l — v^m2) 



(48) 



(49) 



(50) 



These bounds are exact. The optimal structures that realize the bounds are again the wheel assem- 
blages described above. The calculation is similar to those shown above with obvious change in notations. 
Notice that the first term of the bound B2 (|48p corresponds to the bulk stress in the central circle, which 
is proportional to bulk modulus K2 of the material in the circle, and the second term corresponds to the 

stress in the spikes, which is inversely proportional to ki + r;i = — the Young's modulus of the first 

El 

material. This dependence is explained by observing that the radial spikes are under one-directional 
radial load that they transform to a pressure on the central circle. Similar structures realize the bound 
B3 (|49|) . For that case, the first term corresponds to the bulk stress in the center circle, and the second 
- to radial spikes. The second term also depends only on Young's moduli of the materials. The spikes 
in the annulus. Figures [T] [3] transport the hydrostatic stress to the central circle, providing optimal stiff- 
ness. The larger fraction of the stiffest material permits for additional outer annulus that increases the 
stiffness. Figure [2] 



Dual problem The dual elasticity problem is addressed in the same manner. It asks for the minimal 
stiffness of three-material mixture, if one component is rigid. The bounds H47p . (|48[) . H49p are valid with 

a change in notations: k is changed to the plane bulk modulus K , K — — and 77 is changed to shear 
1 

modulus fi, fi — —. 

■n 

The same wheel assemblages are optimal. We observe the dual behavior similar to the conductivity 
problem. The trapezoidal domains in an optimal wheel are filled with the rigid material which transmits 
the applied hydrostatic deformation to the middle circle and material in elastic spikes is under a uniaxial 
deformation in the circumferential direction, the normal (circumferential) stress is constant. 
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